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Abstract 

Using data acquired with the CLEO detector at the CESR e"'“e“ collider at y/s = 3.773 GeV, we 
measure the cross section for the radiative return process e^e“ ^ yJ/iI, J/f) ^ ■, resulting 

;u+/r-) xree(J/b) = 0.3384±0.0058±0.0071 keV, ree(J/b) = 5.68±0.11 ±0.13 keV, 
and = 95.5±2.4±2.4 keV, in which the errors are statistical and systematic, respectively. 

We also determine the ratio Tee[i’{‘IS)\/Tee{J/'f) = 0.45 ± 0.01 ± 0.02. 
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The full and dileptonic widths of a hadronic resonance, Ttot and Tee, describe fundamental 
properties of the strong potential [l|. The value of Tee for a particular resonance is, in 
principle, predictable within QCD, although the strong interaction effects in the quark- 
antiquark pair annihilation make calculations challenmng. Heavy quarkonia offer the best 
testing ground for lattice-based (LQCD) techniques p, and a fortuitous convergence in 
precision to the few percent level is occuring on both the theoretical and experimental 
fronts. In 2003, BaBar measured x Tee 0 , where = B{J/'ip —>• using 

the novel technique of counting radiative returns from e’''e“ collisions at ^/s = 10.58 GeV. 
This allowed the world average 0 uncertainty on TfJJ/'ip) and TtotiJ/i^) to be reduced by 
nearly a factor of two when combined with a BBS 0 determination of the dileptonic 
branching fraction Bn, which has a relative 1.7% uncertainty. At the same time, progress on 
predicting and measuring dielectronic widths of bottomonium 0 0 is occuring, providing 
further checks of LQCD computations. For the T system, the LQCD predictions for the 
ratios ree[T(nS')]/ree[T(lS')] are expected to be more accurate than those of the absolute 
individual widths. 

In this Article we describe a measurement of the J/ijj full and dielectronic widths with 
CLEO-c data from collisions near the peak of the ^’(3770) resonance. The method 

is similar to that used earlier by BaBar 0 and in CLEO’s recent measurement 0 of 
ree[V’( 2 S')]: we select /i+/i“( 7 ) events, each with a dimuon mass in the general region of the 
J/ijj, and count the excess over non-resonant QED production, —>• 7 p+/i“. (The e+e“ 

hnal state is not used due to the large f-channel contribution, which limits the attaintable 
statistical precision relative to The J/ijj component will peak at ^~)=Mji,p 

with a mass resolution dominated by detector effects. The cross section for the excess is 
proportional to B^^ x Tee{J/'ip)- Assuming lepton universality, we can then divide by CLEO’s 
own Bn 0 , with a relative accuracy of 1.18%, once to obtain Tee{J/'il^) and once more for 
rtot(W). 

We use e+e collision data collected with the CLEO detector jof acquired at a center- 
of-mass energy yCo=3.773 GeV at the Cornell Electron Storage Ring (CESR) The 
CLEO detector features a solid angle coverage of 93% for charged and neutral particles. 
The charged particle tracking system operates in a 1.0 T magnetic field along the beam axis 
and achieves a momentum resolution of ~0.6% at momenta of 1 GeV/c. The integrated 
luminosity (£) was measured using 77 , and events [ll| and normalized with 

a Monte Carlo (MC) simulation based on the Babayaga 113 generator combined with 
GEANT-based [l3| detector modeling. Results from the three final states are consistent and 
together yield C = 280.7 ± 2.8 pb“h 

The differential cross section for e+e“ —yJ/V’ —>■ 7 /^^/^” can be expressed 0,1110 
in terms of the e’''e“ invariant mass s, the dimuon mass-squared s', and the variable x = 
(1 — s'/s) as 

^{s,x) = W{s,x) xb{s') xTeeX Bf^f, , (1) 

ax 

where W{s,x) is the initial state radiation (ISR) 7 -emission probability, b{s') is the rela¬ 
tivistic Breit-Wigner function, and Tee is the J/'ip e^e~ partial width (including vacuum 
polarization effects). The ISR kernel, to lowest order in the fine structure constant a, is 

lT(s,x) = — An^ - l^J fl-x +, ( 2 ) 

Tix \ mj J \ 2 J 
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( 3 ) 


in which rUe is the electron mass. The Breit-Wigner function is 6(s') = B{s')/B^^Tee-i 

rj, /N ^ 127rg^^reertot 

(s'-M2)2 + M2rL ’ 

where Ttot is the full width and M the J/'ip mass. 

The cross section cxo = o'(so) for e’''e“ —*• 7J/'0 —»• 7 /r’''/r“ over a specified dimuon mass 
range can be obtained from Eq. (1) and measured: 


(To = 


^J / Ip ^hgd 

ex C 


Tee X ^ 


( 4 ) 


in which is the number of signal events counted, iVbgd is the estimated background, e 
is the detection efficiency obtained from Monte Carlo (MC) simulation, Jq = /(sq), and the 
integral 

I{s) = J W{s,x) b{s') dx (5) 

is effectively insensitive to the value of Ttot- Hence a measurement of ag can be combined 
with measurements to yield For these equations to work, the number of 

events N and the integral above must both be determined with the same limits on x, which 
means the same limits on muon pair mass. 

The above treatment ignores interference effects with the QED production, which 

modify the dimuon mass lineshape asymmetrically around the peak; these effects are in¬ 
cluded at lowest order by replacing B{s') with 


B'{s') 


Aira'^ 

3s' 


QMT tot 

M^-s'- tMT tot 



( 6 ) 


where Q = 3^BeeB^^/a. 

Evaluation of Eq. (6) shows that interference is constructive above the peak and sym¬ 
metrically destructive below, inducing a change in the Breit-Wigner cross section at 
^/s = Mj/^ ± rtot/2, for example, of about ±8%. At ±300 MeV from the peak, the 
interference term induces an effect equal to ±0.8% of the nonresonant QED cross 

section. 

The integral /(s) is performed numerically. The result for our t/s=3.773 GeV dataset and 
a:=0.139-0.488 (i.e., M(/i+/i“)=2.8-3.4 GeV) is /o=188.8±1.3 pb/keV using W{s,x) from 
Eq. (2) and Jo=185.8±1.3 pb/keV when including the radiative corrections in Eq. (28) of 
Ref. Both values for Iq include a net relative increase of ~3% due to the interference 
term (using Eq. (6) instead of Eq. (3)), which occurs because the ISR kernel’s 1/x term 
weights higher masses (constructive interference) more than lower masses (destructive inter¬ 
ference). The quoted uncertainties on Ig are based only on the statistics of the numerical 
integrations. 

The event selection procedure is straightforward. The two highest-momentum 
tracks are required to have opposite charge, individually satisfy either | cos ^1 <0.83 or 
0.85<| cos6'|<0.93 so as to avoid the barrel-to-endcap calorimeter transition region, and 
together have an invariant mass in the range 2.8-3.4 GeV. Bremsstrahlung photons, defined 
as calorimeter showers found within a 100 mrad cone about the initial charged track direc¬ 
tion, are added to the corresponding Lorentz vector for the computation for each 

event. Muon pairs are loosely selected, and electrons effectively vetoed, by requiring the two 
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tracks to satisfy muon-like requirements on the matched energy-to-momentum ratio E/p\ 
the larger of the two E/p values must be <0.5 and the smaller <0.25; electrons typically have 
E/p ~1 and consistently satisfy E/p >0.5. Cosmic rays are suppressed by requiring the pair 
of tracks together to point to within 2 mm of zero in the plane perpendicular to the beams 
and within 40 mm of zero along the beam direction. Cosmic rays are further suppressed by 
requiring the candidate J/^|J to have momentum p^^ = 0.1-1.5 GeV/c, a restriction which 
has no effect upon signal efficiency. 

The dominant backgrounds to a 'yJ/'i/ signal are radiative returns to ‘ip{2S) with subse¬ 
quent decays 'ip{2S) —>• XJ/xjj. Such events will have a true J/i/ but they will 

also tend to have extra tracks or showers, as well as a significant mass recoiling against the 
muon pair. Three requirements are imposed to suppress these events: one on extra charged 
tracks, a second on extra calorimeter showers, and a third on missing mass. The number 
of tracks satisfying loose quality criteria is required to be exactly two. The missing mass, 
|Pcm — +Pm“)I’ where Pcm is the initial state center-of-mass Lorentz vector, and p^+ and 

p^- are the two muon Lorentz vectors, is required to be less than 500 MeV; this value is 
set by the need to reject both charged and neutral titiJ/'i/ events. We search for the most 
energetic shower unassociated with a charged track that is not within a 100 mrad cone of 
either the initial momentum direction of either track or of the opposite of the net muon 
pair momentum direction, and demand it to have energy below 150 MeV. Figure ^ shows 
the two variables sensitive to backgrounds from radiative returns to 'ip{2S) for signal and 
^'ip{2S) —>■ 'yXJ/'ijj, J/'ip —>■ MC, demonstrating that the restrictions on extra showers 


and missing mass separate the signal from these backgrounds. 

The EvtGen event generator jl^, which includes final state radiation (FSR) Hi , and a 
GEANT-based 0 detector simulation are used to study the radiative return (to J /'tp and 
'ijj{2S)) processes with exactly one ISR photon. Events are generated with the polar angle 
distribution from Ref. and account for ISR according to Eqs . (1-3). Non-resonant events 
of the type 7 ( 7 ...)/rare generated from the Babayaga m package, which, unlike the 
radiative return process of Evtgen, includes the effects of multiple ISR and/or FSR photon 
emission, and hence of higher orders in a. The Babayaga code normally includes the 
interference effects with J/ip decays, but this feature was removed for the results shown 
here. 

In the absence of an MG generator package incorporating multiple photon emission from 
the initial state and from J/ip decays in the radiative return process, we calculate the 
efficiency of the J/p) signal in three steps, assuming that ISR and decay radiation are 
factorizable. For the first step, in order to simulate production of 7 ( 7 ...) J/a/ —> 7 ( 7 ...)/i+p“, 
we use a subset of BABAYAGA-generated e+e“ —> events with FSR disabled and 

a muon pair mass (including photons emitted within 100 mrad of either muon’s direction) 
restricted to within ±10 MeV of Mj/^. After detector simulation and reconstruction, 73.2% 
of such events pass our selection criteria. For the second step, we compare the EvtGen 
efficiency with J/p: decay radiation from PHOTOS 3| to that without it, finding that decay 
radiation reduces the efficiency by factor of 0.968. The third step accounts for imperfections 
in modeling track-finding and decay radiation, for which we correct the efficiency by the 
factor 0.995 0, arriving at e=70.5%. 

In order to probe both the background levels as well as the modeling of the largest 
shower and missing mass restrictions, we also perform the analysis with two alternate sets 
of selection criteria, differing from nominal only in that for the “loose” (“tight”) set, we 
require the highest energy shower to have energy less than 200 MeV (100 MeV), and that 
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the missing mass be less than 550 MeV (450 MeV). These variations result in a relative 
efficiency change of +1.6% (—1.9%). 

TablelHshows the expected number of background events from ■j'ip{2S) —>• ■jXJ/'ip, J/'ijj —> 
and from yJ/V', J/V’ — 7 r''' 7 r“. The total number is 1.3% of the signal. For the 
alternate “loose” (“tight”) selection, the relative background prediction is 5% (0.5%). Other 
J/ijj decay modes are found to contaminate our sample at negligibly small levels. No other 
processes will produce a peak at the J/'ip mass. Other backgrounds are assumed to be 
smooth in and httable by a low-order polynomial. 

In order to avoid depending on near-perfect MC simulation of the mass resolution, an 
alternate procedure is used to generate an accurate expected shape of the dimuon mass 
spectrum. We take a clean sample of J/'ip decays from data in which there is 

no interference and convolve the measured mass resolution with the expected effects from 
interference to obtain the expected shapes. CLEO has already accumulated a large sam¬ 
ple of essentially background-free radiative return events from y/s=3.773 GeV to '0(25'), 
'iIj{ 2S) —>■ 7171J/'ll) —>• 7 r 7 r/r+/r“ j^, with almost the same selection criteria as for this analy¬ 
sis. After rejecting events failing the 150 MeV unaffiliated shower veto, we take the mass 
distribution from these 11,305 events, summed over both the charged and neutral dipion 
samples, and offset it by Mj/^ so as to be peaked at zero. The resulting distribution is 
taken to represent the mass resolution function in the 2.8-3.4 GeV mass region. In a toy 
MG, three different mass distributions are generated for dimuons with ISR: from a J/0 de¬ 
cay alone (Eq. (3)), from non-resonant hrst-order QED, e+e“ —> (47ra^/3s'), and for 

the combination including interference ( 47 ra^/ 3 s' + B'{s')). For each “event” in each dis¬ 
tribution, the mass is smeared according to the mass resolution function from y' 0 ( 2 S') data 
and recorded in a histogram. The hnal step is the subtraction of the properly normalized 
QED-only nonresonant mass distribution from that for QED-plus-J/'0-with-interference to 
obtain the expected shape of the J/'ip mass peak from radiative returns. We designate this 
expected shape of the J/ip peak including resolution and interference as H{Mi), which has 
the property H{Mi) = 1 , where Mi is the center of the mass bin. 

The following approach is taken for htting the smooth non-resonant background: in seven 
hts from the widest window ( 2 .8-3.4 GeV) to the narrowest (3.06-3.14 GeV), the lowest 
order polynomial is used that, in combination with the signal shape from above, gives at 
least a 1.0% conhdence level (GL) for the £t. In practice, this meant using a third order 
polynomial for ranges wider than 2.9-3.3 GeV, a linear background for ranges narrower 
than 3.03-3.17 GeV, and a second order polynomial otherwise. We chose this strategy 
to allow for statistical fluctuations in the signal resolution function, to avoid introduction 
of an unphysical background shape, and to maximize the orthogonality of the signal and 
background functions. The htting function is 


/(M,) = Nj/^H{Mi) + AMj2 (7) 

j=0 

where Nj/^ and aj are the hoating £t parameters, AM is the bin width, and, depending on 
the range, 02 and/or 03 can be set to be zero. 

The normalization scheme represented in Eq. (4) will be referred to as the efficiency (E) 
method. The important features of the efficiency method are that one makes no assumptions 
about the composition of the non-resonant background, and that it requires an absolute ef- 
hciency measurement for the selection applied as well as absolute luminosity. An alternate 
normalization scheme, which will be referred to as the ratio (R) method, can be employed 
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to probe several systematic effects. In the ratio method (which was used in Ref. S) , in¬ 
stead of using luminosity measured with 77 , and events, we use the number 

of non-resonant e+e“ — 7 /i^p“ events underneath the J/'ip signal in the muon pair mass 
distribution. This alternate integrated luminosity is 


^ _ ^QED 

—-7-, 

^QED OqeD 


( 8 ) 


in which A^qed is the number of non-resonant y/i+p” events per unit mass from QED alone, 
evaluated at the J/ijj mass (which in terms of our fit is by definition the parameter oq from 
Eq. (7)), eQED is the efficiency for non-resonant QED muon pair events to pass the selections, 
and 5 qed is the cross section per unit mass predicted by non-resonant QED alone, without 
detector effects or selections, at M=Mj/^, including hnal state radiation and vacuum polar¬ 
ization effects. The value for hgEo is obtained by htting the non-resonant —>• 

MC mass distribution at the generator level with the polynomial portion of Eq. (7) only, and 
dividing by the effective luminosity of the MC sample; we find 5qed=0.8510±0.0031 pb/MeV 
for the Babayaga generator and a value 1.7% smaller when using a hrst-order (at most, 
one photon per event) generator [ 2 ^. We compute cqed as in the first and third steps of 
that for e, but with Babayaga FSR enabled. Ending eQED=69.2%. 

The ratio method has, in place of Eq. (4), 




bgd 


Oq 



Tee ^ ^ *^0 ■ 


(9) 


The ratio cqed/ e is expected to be close to unity; systematic eEects which mostly cancel in 
this ratio include those from trigger, reconstruction, radiative corrections, and event selection 
variable modeling. The ratio method replaces systematics associated with the standard 
luminosity measurement and absolute efficiency determination with those applicable to oq 
and (5qed- We note that Nj/^ and oq are almost completely anti-correlated, for which we 
account in the uncertainty propagation. 

Unlike the efficiency method, the ratio method requires understanding non-J/'ip back¬ 
grounds, i.e., the extent to which other final states besides radiative muon pairs populate 
the non-resonant entries in the mass distribution. Based upon measured cross sections 21 


and MC studies of low-multiplicity hadronic final states, we conclude that all other back¬ 
grounds are negligible and assign an uncertainty of 0.3% in x T^e from this source. 

Fits over all mass ranges strongly prefer the shape that includes interference to those 
which do not. This preference for the 2.9-3.3 GeV fit range, as shown with interference in 
Fig. Ha) and without in Fig. Hb), is a 6.6ct effect, as determined from the difference in 
log-likelihoods from the respective fits. The no-interference fit also systematically underes¬ 
timates the yield by ~3%, which, not coincidentally, is the amount by which interference 
changes the overall rate from 2 .8-3.4 GeV. 

Table HH lists the quantities relevant to the i^^^xTee measurements. Gentral values and 
statistical uncertainties of Nj/^ and Oq are taken as unweighted means over the seven fits 
previously described, which have GL’s ranging from 1-18%. Gombining information from fits 
over different mass ranges in this manner samples different relative weightings of background 
and signal regions. The systematic errors on Nj/^ and Oq are taken as the rms spreads of the 
corresponding fit results. The values from the efficiency and ratio methods are consistent 
within their uncorrelated uncertainties. The polynomial fits and shape agree very well with 
the luminosity-normalized expectation from the radiative muon pair MG. 


7 







For the efficiency method, the “loose” (“tight”) selections induce changes in x Fg 


from nominal of +0.5% (—0.0%); similarly, the ratio method variation is +0.0% (+0.3%). 
These small changes demonstrate a good understanding of efficiency and background levels, 
which are reflected in their systematic errors in Table IHl 

The systematic errors on efficiency for the dimuon pair arise by extrapolating errors 
from the CLEO 'ip{2S) XJ/tlj analysis 1^| as appropriate. Quoted errors for J(s) and 
(5 qed include, in addition to MC statistical uncertainties, contributions of 0.5% to account for 
accuracy of the underlying formulae. The systematic error on iVj/p (or Nj/^/ao) from htting, 
computed as described above, is 0.9% (or 1.8%). The accuracy of the htting assumptions 
is tested by pursuing the htting procedure using MC 7'0(2S') —>• 'y7i^7i~fi~ events for the 
htting shape, MC 'yJ/'<p events, without any interference, for signal, and MC for 

non-resonant background, all with statistics much larger than those of the data. From these 
high statistics samples, we hnd that the above procedure introduces a bias in x Fee of 
1 .0% for the efficiency method and 1.1% for the ratio method, both in the upward direction. 
We fully correct for these biases and assign a 0.7% systematic error to the corrections. The 
relative uncertainty attributable to the statistics of the resolution function in the data was 
estimated as follows. An ensemble ■j'ip{2S) MC samples was formed, each of 

the same size as the corresponding sample from the data, and each was used as the resolution 
function for the high-statistics MC muon pair distribution. The relative rms variation of 
the £t results for Nj/^ (oq) for the ensemble was found to be 0.7% (0.8%), which is included 
as an additional systematic uncertainty. If, instead of adding the fitting systematic errors 
described above, we scale up the statistical errors by the square root of the reduced-y^ on 
each of the seven fits, similar total uncertainties are obtained, indicating reasonable error 
assignments. 

Reference asserts that eQEo/e = 1 and that the ratio method is insensitive to radiative 
corrections; we hnd neither to be the case for our analysis due to FSR effects. The ratio 
5 qed/7o changes by 5% when ignoring radiative corrections: when both ISR and FSR are 
allowed, non-resonant events move from higher muon pair mass to lower mass, resulting in 
a net increase in events at M = Mj/^, whereas decay radiation in J/'ifj can only 

shift events out of the signal peak to lower masses. 

We compute an error-weighted average of the E and R methods, accounting for corre¬ 
lations. While the ratio method avoids some of the systematics of the efficiency method, 
it suffers a larger htting error, because A^j/p is almost fully anti-correlated with oq (which 
is the non-resonant muon pair level at This average gives relative weights of 8:1 

for E:R. The weighted average is i3^^xFee=0.3384+0.0058+0.0071 keV. For we use the 
CLEO measurement ^ i3££=(5.953+0.056+0.042)%, yielding Fee=5.68+0.11+0.13 keV, and 
rtot=95.5+2.4 + 2.4 keV. In all cases the hrst errors quoted are statistical and the second 
systematic, and the distinction between the two has been preserved in the propagation of 
uncertainties in Bii^ to Fee and Ftot- 

This measurement of B^^ x Fee is consistent with the BaBar value, and the values 
determined here for Fee and Ftot are more precise and somewhat larger than all previous 
measurements. 

In summary, we have used the radiative return process e^e ~to measure B^^xFee 
with a 2.7% relative uncertainty, and combined this with a CLEO measurement of Bu to 
obtain Fee (3.0%) and F^t (3.6%) with improved precisions. Combining with Fee['0(2S')] = 
2.54+0.03+0.11 keV from Ref. [7| and accounting for common uncertainties of luminos¬ 
ity, Bii^ and lepton tracking, we determine Fee['0(2S')]/Fee(T/'0)=O.45+0.01+0.02 (5.0%), a 



quantity which might be more precisely predictable in LQCD than either Fee alone. 
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TABLE I: Number of background events expected. 


Source 

^ Events 

yJ/'i/), J/ijj 

17.2 

7'0(25) ^ 77r+7r“ J/'0 —> ii~ 

6.2 

7'0(25) ^ yyr'^Tr^ J/'0, 

62.4 

7V’(25) ^ irjJ/'ip, Jji) fi+fj,- 

2.1 

jip{2S) 77r*^J/'0, J/^l^ 

12.0 

7V’(25) ^ 7Xc0, Xco ^ iJ/i’, J/i’ 

0.5 

j'ip{2S) jxci, Xci J/'^ 

18.6 

7V’(25) ^ xXc2, Xc2 iJ/'ik, J/'^ 

49.6 

Total 

168.7 


TABLE 11: Intermediate and final results for the efficiency (E) and ratio (R) methods, with statis¬ 
tical and systematic errors. 


Value 


Quantity 


-^bgd 

Jo (pb/keV) 

Efficiency (E) method 

e(%) 

£ (pb-i) 

Bias Eactor 

^ Tee (keV) 
Ratio (R) method 
ao (MeV-i) 
eQEo/e 

'^QED (pb/MeV) 
(pb“^) 

Bias Eactor 

^ bee (keV) 


12742 ± 202 ± 143 
169 ±9 
185.8 ± 1.6 

70.5 ± 1.2 
280.7 ±2.8 

0.990 ± 0.007 
0.3385 ± 0.0054 ± 0.0075 

165.7 ± 1.6 ±2.1 
0.981 ± 0.008 
0.8510 ± 0.0052 

281.5 ±5.3 
0.989 ± 0.007 

0.3373 ± 0.0087 ± 0.0096 
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FIG. 1: Distributions, for M^^=3.05-3.15 GeV, of the largest calorimeter shower energy unaffiliated 
with a charged track (top) and missing-mass-squared (bottom) for the data (filled circles), in the 
signal yJ/V’ MG (dotted line histogram), 7^/)(25) ^ yVJ/'i/i MG (dashed), and their sum (solid). 
Arrows show the nominal upper limits of values accepted by the event selection. 


11 













3.05 3.10 3.15 

M(|j,+n-) (GeV) 


FIG. 2: Fit (solid line) of the muon pair invariant mass data (filled circles) to the sum of the 
expected shape (dashed line) for a J/'0 decay (a) with and (b) without interference combined with 
a smooth background (second order polynomial). 






